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Nonequilibrium Perturbation Theory 



The goal of this chapter is to construct the perturbation expansion for the 
1-particle contour-ordered Green's function. First, we explain why the usual 
perturbation expansion on the Feynman contour (the real-time axis from — oo 
to oo) or on the Matsubara contour (a segment on the imaginary-time axis 
from — i(3 to i(3) fails in general nonequilibrium situations. 

1.1 Failure of Conventional Time-Ordered Perturbation 
Theory 

The central goal of nonequilibrium many-body theory is to calculate real-time 
correlation functions. For example, we might want to calculate the 1-particle 
time-ordered Green's function, 



tG(x,t;x',f) = (TWx^V,*')]) = TrpT[^(x,^ t (x / ,0] (1-1) 



in the Heisenberg picture, where p is an arbitrary nonequilibrium density 
matrix and the Hamiltonian H(t) is in general time-dependent. 

1.1.1 Equilibrium Many-Body Theory 

In conventional equilibrium many-body theory, we know how to setup a per- 
turbation theory to calculate this quantity. At zero temperature, the density 
matrix is 



where |$o) is the exact ground state of a full interacting but time-independent 
Hamiltonian H. Then the time-ordered Green's function 



(1.2) 



zG(x,t;x',t') = <!PbTO(x,t)^V,0]l#b) 



(1.3) 



is given by conventional Feynman-Dyson perturbation theory, 



2 1 Nonequilibrium Perturbation Theory 



* G (x, t; x',o = (L4) 

(<P |5(oo,-oo)|<P > 

where the caret denotes operators in the interaction picture with respect to 
a quadratic Hamiltonian H where H — H + V, \$o) 1S the ground state of 
Ho, and the S'-matrix is 

^(oo, -oo) = Texp (^-i J dhVih^j (1.5) 

The power series expansion of the S'-matrix and the subsequent use of Wick's 
theorem generates the usual Feynman diagrams, of which the disconnected 
ones cancel against the phase factor in the denominator. 

At finite temperatures, we use the coincidence in functional form of the 
thermal density matrix 

P-PH 

P =— (1-6) 

where Z — Tre~P H is the partition function, and the evolution operator 
U (t) — e~ lHt to setup a perturbation theory in imaginary time. Here again, 
Wick's theorem can be used since each term in the perturbation expansion 
is an average over a noninter acting density matrix po — e~^ H ° /Zq where 
Zo^Tre"^ . 



1.1.2 Failure of Conventional Techniques 

These techniques seem quite powerful, so why can't we apply them to the 
more general problem of Eq. (1.1)? 

Let us consider the Matsubara technique first. Out of equilibrium, there 
is no such thing as a temperature. As a result, in general the density matrix 
is not of exponential form, so there is no way we can use the Matsubara 
trick which consists in a simultaneous expansion of the density matrix and 
the time-evolution operator allowed by the coincidence in functional form of 
these two operators. 

In the Feynman case, once again the density matrix is not a simple pro- 
jector on the ground state as in Eq. (1.2). 

These arguments are correct, but it is instructive to see in more detail 
where exactly does the mathematical construction of the above equilibrium 
perturbation expansions fail. We will try to construct an expansion of the 
nonequilibrium Green's function using the ordinary Feynman approach and 
see that it fails. 

Consider a generic density matrix 



(1.7) 
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where the |<£) can be arbitrary quantum states. First recall that in the Heisen- 
berg picture, the density matrix does not evolve in time since its time evolu- 
tion, given by the quantum Liouville equation, goes in a way opposite to that 
given by the Heisenberg equation of motion, so that the time evolution of p 
cancels out altogether. For convenience, we study one state at a time, namely 
we want to calculate the expectation value 

iG#(x,t;x',t / ) = (^T^Cx^V, (1-8) 

and we can obtain the full Green's function by G — ^2<pP<pG&. 

Now consider the following partition of the Hamiltonian H(t) — H + H'{i) 
where H is the unperturbed equilibrium Hamiltonian (but may still contain 
interactions) and all the time dependence is included in the nonequilibrium 
perturbation H'(t). The field operator in the interaction picture is 

4>(pc,t) = e im rl>(pt)e- iHt (1.9) 
and in the Heisenberg picture, 

V>(x, t) = U^{t)ip(x)U(t) = 5(0, £)V>(x, t)S(t, 0) (1.10) 
where S(t, 0) = e lHt U(t) and the evolution operator is 

U(t) =Texp(-i dhnih)^ (1.11) 

more generally, we have the /S-matrix 

S{t,t') = e iHt U(t,t')e- iHt ' = Texp ^-i jf dtxH'it^ (1.12) 

where U(t,t') = Texp (-i f* dtiH(h)j and H'(t) = e iHt H' \t)e~ iHt is the 

nonequilibrium perturbation in the interaction picture 1 . 
We substitute these relations into Eq. (1.8), 

iG$(x,t;x. f ,t f ) = (^\T[S(0,t)^t)S(t,0)S(0,t')^(^,t f )S(t',0)]\^) (1.13) 

Now the interaction picture state is given by 

\$(t))i = S(t,0)\*(P)h =S(tM*) (1-14) 
since all pictures coincide at t = 0. We can then say that 

1 We use curly letters for U and S to indicate that these operators take care of the 
full nonequilibrium Hamiltonian TL with H' as the perturbation. Later upright 
letters U and S will be used for the corresponding operators taking care of the 
equilibrium Hamiltonian H with whatever interactions it may contain as the 
perturbation. 
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|<£) = 5(0, ±oo)|<2>(±oo))/ (1.15) 
Putting everything together, we have 

iG*(x,£;x',t') = (<P(oo)|/5(oo,0) 
xT[<S(0, £)^(x, t)S(t, t / )^ t (x / , *')£(*', 0)]5(0, -oo)|<P(-oo))j (1.16) 

Now, the S'-matrix <S(oo, — oo) is itself a time-ordered product, so that we can 
write 

iG*(x,t;x',t') = (^(oo)| / T[<S(oo,-oo)^(x,t)^ t (x / ,t / )]l^(-oo))/ (1-17) 

making use of the group property of S(t, t'). 

So far everything is correct in a general nonequilibrium setting. At this 
point, the conventional Feynman expression Eq. (1.4) is obtained only provided 
that the following crucial step holds: |<P(oo))j and |<P(— oo))j differ only by a 
phase factor. This is true in equilibrium but breaks down out of equilibrium. 

The reason is the following. In the conventional zero-temperature theory, 
we adiabatically switch on and off the interaction: 

H e = H + e~ eW H' (1.18) 

where e is a positive infinitesimal. We also choose |<£) to be |$b), the ground 
state of H — H + H' . We assume this ground state to be nondegenerate. Since 
the time evolution is adiabatic, we can use the adiabatic theorem to say that 
|<P(±oo))j = lim e ^o <S e (±oo, 0)|$b) are both eigenstates of H. Since |#b) is 
nondegenerate, these two states can differ by at most a phase factor, 

|<2>(oo))j = e iL |#(-cx)))/ (1.19) 

For definiteness, using the notation |<£q) = |^( — °°))/ it is eas Y to see that 

e lL = <<2> |<S(oc,-oc)|<Z> ) (1-20) 

so that Eq. (1.17) gives Eq. (1.4) directly. 

Out of equilibrium, the use of the adiabatic theorem is unjustified since 
under the assumption of a general time-dependent Hamiltonian H(t) the 
time evolution is not adiabatic. In the case that the Hamiltonian is time- 
independent but the density matrix is still arbitrary - the important case of 
a nonequilibrium steady state, we still cannot use the adiabatic theorem be- 
cause the generic quantum states \<1>) are not necessarily eigenstates of H and 
the adiabatic theorem applies only to eigenstates of the Hamiltonian. If we 
insist and expand the density matrix in the basis of eigenstates of H, it is in 
general not diagonal in this basis, 

Tr^O = Y^P^PW) = $>„„'<n|0|n') (1.21) 

<P nn' 
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so that even if |n(oo))j = e lLn \n{— oo))j, the appearance of crossed terms 
spoils the expansion. Finally, even if the ground state is assumed to be nonde- 
generate, excited states \n) appearing in the nonequilibrium density matrix p 
can be degenerate so they can be mixed by a non-Abelian Berry phase under 
adiabatic evolution, which once again invalidates the conventional procedure. 

In brief, out of equilibrium |<£(oo))j and |<P(— oo))j are not simply related. 
This is obvious physically since when a system is driven out of equilibrium, 
its asymptotic future is quite different from its initial preparation in the re- 
mote past. The purpose of the various nonequilibrium contours such as the 
two-branch Schwinger-Keldysh contour is precisely to avoid any reference to 
the state in the asymptotic future |<£(oo))j and base the expansion solely on 
\$(—oo))i, the state in the asymptotic past. 

1.2 The Schwinger-Keldysh Contour 

This idea was pioneering by Schwinger [1], Craig [2], Mills [3], and popularized 
by Keldysh [4]. To get rid of the unwanted asymptotic future state \<P(oo))i 
we simply use 

|#(oo))j = S(oo,-oo)|#(-oo))/ (1-22) 

namely, we rewind the time evolution back to the asymptotic past. Equation 
(1.17) then becomes 

zO^x, 1; x', if) = (<2>(-oo)| jS(-oo, oo)T[<S(oo, -oo)^(x, £)^ f (x', if)] |#(-oo)) / 

(1.23) 



£^ 

C- 



C = C+UC- 



Fig. 1.1. Schwinger-Keldysh contour C. 
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Now we cannot simply push <S(— 00,00) past the time-ordering operator T 
and merge it with the forward evolution S (00, —00), since the whole backward 
evolution S{— 00, 00) lies to the left of the T-product. Left means later as far 
as time ordering is concerned. A clever way to kill these two birds with one 
stone is to introduce a two-branch contour C and ordering along this contour 
enforced by a contour- ordering operator T c . This Schwinger-Keldysh contour 
(Fig. 1.1) consists of two oriented branches C — C + UC_, the forward branch 
C+ extending from —00 to 00 and the backward branch C_ extending from 
00 to —00. We now extend the time variables t, t' to variables defined on this 
contour r, r', and define the fundamental object of nonequilibrium many-body 
theory, the contour-ordered Green's function, 

zG(x,t;x',t') = <T c [^(x,TtyV,T')]> = Tv pT c [^r)^ ^.\r')\ (1.24) 

With variables and ordering along the contour C, the interaction picture ex- 
pression corresponding to Eq. (1.23) is 

iG*(x,r;x / ,r / ) = <<P(-oo)|iT c [S c (-oo, -oo)^(x,t)^( x / ,t , )]|^(-oo))j 

(1.25) 

where the contour S- matrix is 




(1.26) 



where the integral is a line integral along the contour C . Note that there is 
no phase factor in the denominator as in the equilibrium theories since the 
contour ^-matrix by itself (i.e. not in a T c -ordered product) is actually unity. 
According to our previous discussion, the full Green's function is now given 

by 

iG(x, r; x', r') = Tr p(-oo)T c [S c (-oo, -oo)^(x, t)^(x', t')] (1.27) 
where we have a density matrix defined by 

p(-oo) = £p*|#(-oo))/<#(-°°)l/ (1-28) 

where now the states are not the Heisenberg states |<£) but the interaction 
picture states in the remote past |<£(— oo))j. From Eqs. (1.7) and (1.15), it is 
not hard to show that 

p = S(O,-oo),0(-oo)<S(-oo,O) (1.29) 

So far this is rather general. Now, in the Keldysh theory we do adiabatically 
switch on the nonequilibrium perturbation V according to e et from the remote 
past t — —00 to the present t = 0, but we do not switch it off. This enables us 
to study stationary nonequilibrium states, a long time after the system has 
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been driven out of equilibrium 2 . Because of the adiabatic switch-on procedure, 
p{— oo) can be identified as the density matrix of the system in the remote past 
when the nonequilibrium perturbation is turned off. It is then evolved in the 
usual way by the S'-matrix <S(0, — oo) to a nonequilibrium density matrix p. 
One thus usually chooses p{— oo) to be an equilibrium distribution, p{—oo) — 



1.3 The Closed Time Path Contour 

In the previous section, we constructed the perturbation expansion such that 
the Green's function would be expressed as an average over the density matrix 
of the system in the remote past because this is what we know: because of the 
adiabatic switch-on procedure, we have p(— oo) = e _/3H jZ which is explicitly 
known. Indeed, the idea of the Keldysh technique is that we want to avoid 
reference to the general nonequilibrium density matrix p which is not usually 
explicitly known. However, this is not mandatory. If we actually know the 
density matrix p of the system at t — 0, then the early Eq. (1.13) is already 
enough: 

iG(x, r; x', t') = Tr pT c [S c (0, 0)^(x, r)^V, t')} (1.30) 

where p is now the full density matrix Eq. (1.7) at t — 0, and the contour over 
which the T c operator and the S'-matrix S c are defined is the so-called closed 



Co 



max(t,t ; ) 



to 



min(t, t') 



Fig. 1.2. Closed time path contour Co- 

2 Actually, we may still study transient behavior in the Keldysh formalism, but 
provided we use an explicitly time-dependent Hamiltonian and there are no initial 
correlations, namely, the initial density matrix has to be noninteracting. 
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time path contour Co going from to the latest of t or t' , and then back to 
0. Furthermore, there is nothing special about the time t — 0: we might have 
defined the Heisenberg and interaction pictures with respect to some other 
initial time to. Then the average would be over the density matrix at time to: 



and the closed time path contour [5, 6] now goes from to to the latest of 
t or t', and then back to to (Fig. 1.2). To recover the Schwinger-Keldysh 
contour, we take the limit to — > — oo and insert a factor of S(t, oo)<S(oo, t) 
or S(t', oo)<S(oo, t') in the perturbation expansion Eq. (1.13) to extend the 
contour to oo, depending on whether t or t' is the latest. The density matrix 
is now p{— oo) and we invoke the adiabatic switch-on procedure to choose 
p(-oo) = e -P H /Z. 

1.4 Interactions: the Kadanoff-Baym Contour 

Equation (1.31) is in a form suitable for a perturbation expansion provided 
that p(to) is a noninteracting (i.e., 1-particle) density matrix and that the 
field operators if), ^ in the interaction picture evolve with a noninteracting 
Hamiltonian H. Indeed, these are the conditions of applicability of Wick's the- 
orem [7]. They are satisfied in the Keldysh theory if H is noninteracting since 
then p{— oo) = e~^ H jZ is a 1-particle density matrix. What if H contains 
interactions? 

Let us first keep a general to and assume that p{to) — e~^ H jZ with H a 
general interacting Hamiltonian. We can always take to — > — oo at the end to 
recover the Keldysh theory. We first break further the equilibrium Hamiltonian 
H — Ho + V into a noninteracting part Ho and the interactions V. We want 
to express the interacting density matrix e~ !3H in terms of a noninteracting 
one e _/3ff °. We make use of Eqs. (1.11) and (1.12), but in imaginary time 
and with a general initial time to- Now the full (equilibrium) Hamiltonian 
is H. As mentioned earlier, we will therefore use upright letters U and S 
for the evolution and ^-matrix respectively. The interaction and Heisenberg 
pictures are now defined with respect to the initial time to, meaning that 



A(t) = e iH(t-t ) Ae -iH(t-t ) and A(t) = e iH (t-t ) Ae -iH (t-t ) m g ince jj ig 



time-independent, the evolution operator is simply U(t,to) — e tH( - t *°) and 



iG(x, r; x', r') = Tr p(t )T c [S c (t , t )$ (x, r)^( x ', /)] 



(1.31) 



we have 




. We then have 



(1.32) 



(1.33) 



where here the S'-matrix S(to — ifi, to) evolves the density matrix along a 
contour [to, to — i/3] on the imaginary axis. The Green's function Eq. (1.31) 
now assumes the form 



1.4 Interactions: the Kadanoff-Baym Contour 



9 



iG(l, T) = i Tr e-P Ho S(t - i/3, t )T c [S c (t , t Q )^ H {l)^ H {l')\ (1.34) 

where for simplicity we denote space and time arguments by a collective in- 
dex 1 = (x,r), 1' = (x.',r'). We have also explicitly indicated that the field 
operators still evolve according to the full H, see Eq. (1.9). They have to be 
brought to the current interaction picture defined with respect to Ho. This 
can be achieved as previously through the S'-matrix, 

^h(x, t) = S(f , (x, t)S(t, t ) (1.35) 

where the caret without the subscript denotes the interaction picture with 
respect to H . Equation (1.34) thus becomes 

iG{l,l') = ^TTe-P H °S{t -i(3,to) 

xT c [S c (t ,t )S(t ,t)i>(l)S(t, t )S(t , t')^{l')S{t', to)} (1.36) 

where we now see that the S'-matrix for interactions S evolves along a three- 
branch contour C* — Co U [to, to — i(3] which is the Kadanoff-Baym contour [8] 
(Fig. 1.3). This contour goes from to to the latest of t and t' , back to to, and 



C max(t,t r ) 

to - \ • > 



min(£, t') 

to - ifi 

C* = C U[t 0i t -i(3] 



Fig. 1.3. Kadanoff-Baym contour C* . 



then down to to — ifi- We define an ordering operator T c * along this contour so 
that we can move the thermal ^-matrix S(to — if3, to) into the contour-ordered 
product. We then obtain 

. rn V] = Tr e-^Tc* [S c * (t - tg, t )S c {t , t )ij>(l)ft (!')] 

l ' ) Tre-P H o Tc *[S c *(t -i(3,t )Sc(toM} 1 } 
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where in analogy with Eq. (1.26) we define a new contour-ordered ^-matrix 
defined along the Kadanoff-Baym contour, 

S c .(t -*/Mo) = T c .exp(-*jr dnVinYj (1.38) 
We have written the partition function as 

Z = Tr e~ fiH = Tr e~ 0H ° S{t - %P, to) = Tr e~ 0H °T c * [S c * (t - i(3, t )S c (to,t )} 

(1.39) 

since the S'-matrices are already time-ordered on their respective contours, 
and since times on the [to, to — i(3\ strip are always later than times on Co, we 
simply have <S c (to,to) — 1 and S c *(to — i[3,to) — S(to — i(3,to) even inside the 
T c * -ordered product. Note that the time arguments in Eq. (1.37) are defined 
on the three-branch Kadanoff-Baym contour C* . Equation (1.37) is now fit 
for perturbation theory since by assumption the averages are with respect to 
a 1-particle density matrix e _/3ff ° and the field operators evolve according to 
the noninter acting Ho, so Wick's theorem can be applied. 

The Kadanoff-Baym formalism is adequate for the study of initial corre- 
lations, namely the effect for times t > to of having an interacting density 
matrix at time to, without the assumption t ^> to. The price to pay for this 
powerful formalism is that the Green's function is defined on a three-branch 
contour and has a complicated expression in terms of the simultaneous per- 
turbation expansion of two ^-matrices. However for many practical purposes 
this is overkill, and for steady-state problems we do not care about the effect 
of initial correlations. In many cases we can assume that correlations decay in 
time so that if we take the limit to — > — oo, at any finite time t ^> to there is 
no signature left of the correlations in the initial density matrix p(to). This is 
the Bogoliubov principle of weakening correlations, a general principle in non- 
equilibrium statistical mechanics. It is however advised to keep in mind that 
in some cases initial correlations can persist at long times due for example to 
the presence of metastable states. 

1.4.1 Neglect of Initial Correlations and Schwinger-Keldysh Limit 

For most practical purposes we can safely ignore the initial correlations when 
taking the limit to — >■ — oo. It has been shown [9, 10, 11] that neglecting 
initial correlations amounts to neglecting the imaginary strip [to, to—ifi] in the 
Kadanoff-Baym contour, so we are back with the Schwinger-Keldysh contour 
(after having extended the closed time path contour as explained in section 
1.3). In this limit, the S'-matrices in the denominator of Eq. (1.37) are trivial 
and the denominator is just Tre _/3H °. We thus have 

tG(l, 1') = TipoT c [S c {-oo, -oc)<S c (-oc, -00)^(1)^(1')] (1-40) 

where the noninter acting density matrix is just 
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e -pH 

p0 = T^m; (L41) 

and S c (— oo, — oo) is just the Kadanoff-Baym contour-ordered ,S-matrix of 
Eq. (1.38) but neglecting the third branch of the contour and taking the limit 

t —> -oo, 

5 c (-oo, -oo) = T c exp ^-i j> dr x V(ti)^J (1.42) 

so that now all the ordering takes place along the Schwinger-Keldysh contour 
C . If the system is initially at zero temperature, we simply have 

iG(l, 1') = <<?>o|T c [S c (-oo, -oo)<S c (-oo, -oo)^(l)^(l')]|*o) (1-43) 

with |^o) the noninteracting ground state of Hq. We see that we have now 
almost recovered Eq. (1.27), albeit with an additional S- matrix S c contain- 
ing the effect of the interactions V . We summarize both zero and finite- 
temperature results Eqs. (1.43) and (1.40) in the formula 

iG(l, 1') = <T c [S c (-oo, -oc)<S c (-oo, -oo)^(l)^(i')]) (1.44) 

with the suitable expectation value. This is the essential starting point for cal- 
culations in Keldysh theory. The perturbation expansion can now be carried 
out, with both nonequilibrium terms H' from S c (Eq. (1.26)) and interaction 
terms V from S c (Eq. (1.38)) appearing in the expansion and in the ensuing 
Feynman diagrams. Let us remark here that the power of the Kadanoff-Baym 
and Keldysh approaches to nonequilibrium many-body theory - or more gen- 
erally, to nonequilibrium field theory - lies in its structure being formally 
identical to that of usual equilibrium many-body theory, albeit with a time 
evolution and the corresponding perturbative expansion defined on a more 
general contour. Then most of the tools of quantum field theory can be ap- 
plied: Feynman diagrams, integral equations for vertex functions such as the 
Dyson equation, etc. 



1.5 Contour Dyson Equation 

The whole perturbation expansion for the contour-ordered 1-particle Green's 
function on the nonequilibrium contour C can be resummed in the form of an 
integral equation, the Dyson equation: 

G(1,T) = G (1,1') + J d2G (l,2)U(2)G(2,l f ) 

+ J d2 J d3G (l, 2)17(2, 3)G(3, T) (1.45) 

G(1,T) = Go(l,l') + J d2G(l,2)/7(2)G'o(2,l / ) 

+ J d2 J d3G(l, 2)17(2, 3)G (3,T) (1.46) 
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where G(l,l') = -i(T c [ip(l)^(l')]) is the exact Green's function Eq. (1.24) 
and Go(l, 1') = —i(T c [ip(l)ip^(l')]) is the unperturbed Green's function with 
field operators in the interaction picture, U(2) is a 1-particle potential and 
X"(2, 3) is the 1-particle irreducible self-energy. The integral sign means a sum 
over all internal variables, J d2 = ^ / d*2 J c olr^. For simplicity, we will use 
the following notation, 

G = G + GqUG + GqSG (1.47) 
G = Go + GUG + GUGo (1.48) 

1.6 Initial Correlations with Arbitrary Initial Density 
Matrix 

For the sake of completeness, let us just mention a few words about the 
most general problem (as far as I know) although we won't actually try to 
solve it in these lectures. Consider Eq. (1.31). In the Keldysh theory we let 
to — > — oo and neglected initial correlations. In the Kadanoff-Baym theory we 
considered initial correlations by keeping a finite to, but we assumed that the 
initial density matrix p(to) had the equilibrium form p(to) — e~ f3H jZ with H 
an interacting Hamiltonian. Now, what if p(to) is a general, nonequilibrium, 
interacting density matrix? This corresponds to preparing the system in a 
correlated nonequilibrium state at a given time to and observing the time 
evolution of the system for finite times t > to without assuming t 3> to. People 
need that in the study of correlated plasmas for example. This problem has 
been studied by Fujita [9], Hall [10], Kukharenko and Tikhodeev [11], and 
Wagner [12] who gives a rather clear and comprehensive discussion. The idea 
is to use a modified Kadanoff-Baym contour. First, the general nonequilibrium 
density matrix can still be written in the form 

-AS 

since it is positive definite and Hermitian, where A is not the temperature 
and B is not the Hamiltonian but some general quantities. But because of the 
formal analogy to (3 and H, we can setup a perturbation expansion on a mod- 
ified Kadanoff-Baym contour with the imaginary strip [to — iX, to]. This very 
general approach encompasses the Feynman, Matsubara, Schwinger-Keldysh, 
and Kadanoff-Baym approaches as special cases [12]. 

1.7 Relation to Real-Time Green's Functions 

In all that follows, we will confine ourselves exclusively to the Keldysh ap- 
proach [4, 14, 16, 13, 15] (but the equilibrium Hamiltonian H can still contain 
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interactions). Let us recap. Our initial goal is to calculate the real-time Green's 
function of Eq. (1.1), 

iG(x,*;x',0 = (IWx^V,*')]) =TrpT[^(x,t)^ t (x / ,0] (1-50) 

with operators in the Heisenberg picture with respect to the full nonequilib- 
rium Hamiltonian TL and a nonequilibrium density matrix p. We have defined 
a contour-ordered Green's function Eq. (1.24) with time arguments on a con- 
tour G, 

zG(x,t;x',t') = (TcMx.rWVy)]) (1-51) 
and obtained a perturbation expansion for that Green's function in Eq. (1.44), 

zG(x,t;x',t') = (T c [S c (-oo, -oo)<S c (-oo, -oo)^(x, r)^V, r')}) (1.52) 

with operators in the interaction picture with respect to the noninteract- 
ing part Hq of the equilibrium Hamiltonian H. Now, how is the contour- 
ordered Green's function Eq. (1.51) related to the real-time Green's function 
Eq. (1.50)? 

The answer is that depending on which branch the contour arguments r, r' 
belong to, different real-time Green's functions are obtained: 



G(x,t;x',t') = < 



f G T (x,t;x / ,f) = -i(r[^(x,t)^t( x ' } f)]) if r>T ' G C+, 
G<(x,t;x',t') = i(^t( x ' } f)^( x> t)) if r G C+,r / G C_, 
G > (x,t;x',t / ) = -^(x,*)^*',*')) if t G G_,r' G C+, 
k G f (x,£;x',£') = -i(T[^(x,t)^ t (x , ,t , )]> if t,t' G C_, 

(1.53) 

where G T is the time-ordered Green's function of Eq. (1.50) and T is an anti- 
time ordering operator which orders operators in the opposite way as T, so 
that G T is the anti-time-ordered Green's function. G <J> are the lesser and 
greater Green's functions, respectively. One then defines a 2 x 2 matrix Green's 
function in real time, 

G =|JS) = I (i.54) 

It is not hard to show from the definitions that the following identities hold 



Q R 


— G T 


-G< 


= G> 


— G T 


(1.55) 


G A 


— G T 


-G> 


= G< 


-G f 


(1.56) 


G K 


= G> 


+ G< 


= G T 


+ G f 


(1.57) 



1.7.1 Larkin-Ovchinnikov Representation 



In the so-called Larkin-Ovchinnikov representation [16], we perform a linear 
transformation on G to obtain another matrix G, 
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G = Lt 3 GL1=( G *^ (1.58) 

where L — ^(t° — it 2 ) where t % are the Pauli matrices (r° is unity). G R,A 
are the retarded and advanced Green's functions, 

G R {*,t;x! ,t') = -iO(t-lf){{ip(x,t),ipt(x!,lf)}) (1.59) 
G A (x, t; x', t') = id(t' - t){{tp(pc, t), ^(x', *')}) (1-60) 

and G K is the Keldysh Green's function, 

G K (x,t;x!,lf) = -zMx,t),^(x'^)]} (1-61) 

1.7.2 Langreth Theorem of Analytic Continuation 

Consider the 'matrix products' occurring in Eqs. (1.45) and (1.46). These are 
convolution integrals on the two-branch contour. Consider first the following 
product of two contour-ordered quantities (such as Green's functions or self- 
energies) , 

C(l,l') = J d2A(l,2)B(2,l') (1.62) 
or in simple matrix notation, 

C = AB (1.63) 

How do we convert this product to an integral over the real axis involving 
the real-time components of A and B? This is accomplished by the Langreth 



theorem [18, 19] which consists in a series of rules: 

(AB)> = A R B> + A>B A (1.64) 

(AB) R,A = A R,A B R,A (1.65) 

{ABC)> = A R B R C> + A R B>C A + A>B A C A (1.66) 

{ABC) R > A = A R > A B R > A C R > A (1.67) 



where the 'matrix products' on the right-hand side consist in summations 
over internal degrees of freedom (space and spin) and convolution integrals 
on the real axis from — oo to oo. We are now completely rid of the contour; 
the 2x2 matrix structure takes care of the nonequilibrium physics and we 
can use ordinary real-time integrals. In particular, we can Fourier transform 
to frequency space in the case of a nonequilibrium steady state in which the 
Green's functions are time-translationally invariant, Gij(t,t') — G^{t — t'). 

We will not derive all these rules but only show how the first one can be 
obtained. Consider C — AB on the Schwinger-Keldysh contour C. Then 
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Jc 

/OO p — OO 

dtiA T (t,t 1 )B < (t 1 ,t')-\- / dt 1 A < (t,t 1 )B T (ti,t') 
-oo J oo 

/oo 
dti 
-oo 

From Eqs. (1.55) and (1.56), we know that A T = A R +A< and B f = B<-B A , 
hence 

/oo 
dh [(A R (t,t 1 )+A<(t,t 1 )) B<(t u t') - A<(t,h) (B<(t 1 ,t') - B A (t u t'))} 
-oo 

/oo 
dh [A R (t 7 t 1 )B<(t l7 t') + A<(t, h)B A (t l7 t')] 
-oo 

The other rules are obtained in a similar fashion. Usually, one writes down 
the diagrams as in conventional many-body theory and the corresponding 
analytic expressions involving contour-ordered Green's functions. Typically, 
these involve integrals over the contour. Then one translates these expressions 
in the real-time language using the Langreth rules. This is often referred to 
as 'analytic continuation' from the contour to the real-time axis in analogy 
with the Matsubara formalism. We feel this is somewhat pedantic since the 
Schwinger-Keldysh contour is not really a contour in the complex plane; it is 
only two counter- propagating copies of the real axis. 



A T (t, t 1 )B<{t 1 ,t') - A<(t, t 1 )B T {t 1 ,t') 



2 



Quantum Kinetic Equations 



We have now set up the necessary machinery to construct quantum kinetic 
equations which govern the evolution of correlation functions in a quantum 
many-body system driven out of equilibrium. We first say 'correlation func- 
tions' which are two-time objects instead of 'distribution functions' which 
are one-time objects because equations typically do not close for distribution 
functions. However equations for correlation functions are very complicated 
integral equations which can usually only be solved numerically. There how- 
ever exists a variety of useful approximations which make the equations for 
distribution functions close. 

We will thus first discuss the Keldysh and Kadanoff-Baym equations which 
are exact quantum kinetic equations for the nonequilibrium correlation func- 
tions of the system. Then we will introduce the Wigner transform, the gradient 
expansion and the gradient approximation which yield the quantum Boltz- 
mann equation (QBE). As this equation is still rather complicated, further 
approximations can be made such as the quasiparticle and quasiclassical ap- 
proximations. The QBE is still an equation for a four-parameter (two-point) 
correlation function G < (p, uj, R, T). Expressed in terms of the Wigner dis- 
tribution function f w (p,H, T) = — iG < (p,t — 0,R, T) which is a one-time 
function, the QBE does not close. One can however generate a closed equation 
for the Wigner distribution function by introducing the so-called generalized 
Kadanoff-Baym ansatz [17]. We will also recover the classical Boltzmann equa- 
tion as a limiting case of the QBE. 

We can actually understand from simple arguments why a kinetic equation 
for a quantum distribution function should have more terms than that for a 
classical distribution function. Recall how the classical Boltzmann equation 
is derived. Because of the incompressibility of phase space, the total rate of 
change of the classical distribution function /(p, R, T) is equal to the rate of 

change of / due to collisions, say ( ) = I[f]' 




dT 



I[f] 
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The Boltzmann equation is then obtained simply by using the chain rule for 
derivatives, 

d 



dT 



+ v-V R + F-V p }f = I[f] 



(2.1) 



since dH/dT — v is the velocity and dp/dT = F is the force from Newton's 
second law. In the quantum mechanical case, the distribution function also 
has an energy argument lv: f — f(p, uj, R, T). We thus generate an additional 
term in the corresponding Boltzmann equation: 



d_ 
dT 



+ v- V R + F- V 



du> d 



p + dT J 



f = nn 



But the rate of change of energy is just the power, duj/dT — P — F • v. 
Considering electric and magnetic fields, the force is just the Lorentz force 
F = e(E + v x B), so du/dT — eE • v since the magnetic force does no work. 
Hence we obtain [13] 



d_ 
df 



+ v • ( Vr + eE-|- ] + HE + v x B) • Vj 



/ = I[f] (2-2) 



We will see later in this chapter that this is essentially the structure of the 
QBE. 



2.1 Keldysh Equation 

We first apply the fourth Langreth rule Eq. (1.67) to the contour Dyson 
equations (1.47) and (1.48). We obtain 

G R,A = G R,A + qR,Ajj qR,A + qR,A jjR,AqR,A ^ 
G R,A = G R,A + qRA UG RA + G R,A S R,A G R,A ^ 

because a one-body potential U is proportional to a delta function in time (in- 
stantaneous) and is thus neither retarded nor advanced. Hence we see that the 
retarded and advanced nonequilibrium Green's functions obey Dyson equa- 
tions which are formally identical to the equilibrium Dyson equations. If we 
now apply the third Langreth rule Eq. (1.66), we obtain 

G> = + G§UG> + UG a + G§S R G> + G^S>G A + S a G a 
G> = g| + G R UG% + G>UG A + G r U r Gq + G R S>G A + G>E A G A 

because the instantaneous one-body potential is diagonal in Keldysh space 

< 

(U > — 0). These two equations are the Keldysh equations. Integration over 
the real axis is understood, for example 
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/OO p 
dh / dx 1 G^(x,t-x 1 ,t 1 )U(x 1 ,t 1 )G < (x 1 ,t 1 -x',t') 
-oo J 

and 

/OO /• />oo p 

dti I dxi / cft 2 / dx 2 Gj(x,t;xi,ti) 
-oo </ J — oo J 

x r fl (xi , * i ; x 2 , * 2 ) C< (x 2 , 1 2 ; x' , t') 

The Keldysh equations can be iterated. The infinite order iterate of either 
equation gives the same following explicit equation for G K , 

G> = [1 + G R (U + S R )}G> [1 + (U + U A )G A ] + G R U>G A (2.5) 

which is the usual general form of the Keldysh equation. 

We now digress a moment to study the free Green's functions Gq ,>,r,a 
more closely. They are propagators for a noninter acting Hamiltonian Hq. To be 
fully general, we consider a multiband one-particle Hamiltonian Hq V (x, —iV) 
where fi, v are band (and/or spin) indices. We however require that each Hq v 
be a properly symmetrized function of the conjugate variables, 

i^x, -iV) = ]T {/r(x), ^"HV)} + «^(x) + v^(-iV) (2.6) 

i 

Then Hermiticity of the second quantized Hamiltonian, 

tf = W dx^t(x)^(x, -zV)^(x), (2.7) 

requires that 1 

^(x,-iV) = ^(x,iV)* (2.8) 

We use a matrix notation where Gq V is represented by a matrix Gq in band 
space. It is easy to show by application of the equation of motion tech- 
nique that the contour-ordered propagator Go(x, r; x', t') satisfies the fol- 
lowing equations of motion 2 where matrix multiplication in band space is 
understood, 



1 Two remarks. First, in the case that the only imaginary quantities in Hq v come 
from — iV (for example for a spin system without Zeeman term and in the absence 
of spin-orbit coupling), Eq. (2.8) requires the matrix Hq" to be symmetric. Sec- 
ond, if Hq v depends only on x and not on — iV (for example for a tight-binding 
model where x becomes a discrete variable and V is replaced by finite differences), 
Eq. (2.8) requires the matrix Hq v to be Hermitian. 

2 The delta function 6(r, t') is defined on the contour, such that S R ' A (t, t') = S(t—t') 

< 

and 5>(t,t') = 0. 
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[id T - # (x, -*V)) G (x, r; x', r') = <5(x - x / )<5(r, r') (2.9) 

G (x, r; x', r') (-*0 T , - tf (x', iV')) = <5(x - x')5(t, r') (2.10) 

where in Eq. (2.10) the derivatives d T > and V act to the left. We do not write 
explicitly unit matrices in band space. If we take the real-time components of 
these contour equations, we obtain 

(id t - # (x, -iV)) Cj' A (x, £; x', = 5(x - x')($(* - *') (2.11) 

(z<9, -^o(x,-zV))c o < ' > (x,t;x / ,t / ) = (2.12) 

Gj' A (x, £; x', If) (-id t > - H (x!, iV')) = <5(x - x')<$(* - *') (2.13) 

Gf' > (x,t;x , ,t , )(-^ -fTo(x',iV')) =0 (2.14) 

We can write in simplified notation 

G^G R > A = 1 (2.15) 

Gq 1 Gq = (2.16) 
Gq^Gq 1 = 1 (2.17) 

G> G" 1 = (2.18) 

where the action of the operator Gq 1 on the right and on the left is defined 
in Eqs. (2.11)-(2.14). These relations will be helpful in deriving the Kadanoff- 
Baym and quantum Boltzmann equations in the next section. 

The Keldysh equations are used directly in mesoscopic transport for ex- 
ample, where one studies finite-sized nonequilibrium systems with important 
spatial and/or temporal inhomogeneities. Then the continuous spatial argu- 
ments x, x' in the Green's functions are typically replaced by discrete site 
indices i,j (i.e. in a tight-binding representation) or by principal quantum 
numbers n,n' for discrete energy levels (i.e. in quantum dot physics). The 
Green's functions then become finite-sized matrices and the Keldysh equa- 
tions become matrix equations. For a finite-sized system, it is possible to 
show that the Keldysh equation (2.5) becomes 

G> = G r Gq 1 G^ [1 + (U + S A )G A ] + G R E>G A (2.19) 

but the first term vanishes by Eq. (2.16), so that we have 

< < 

G> = G R S>G A (2.20) 
which is the form of the Keldysh equation that is used in mesoscopic transport. 
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At the other side of the spectrum, if one wishes to study transport in 
macroscopic samples, i.e. transport in metals or bulk semiconductors, non- 
equilibrium (but not mesoscopic) superconductivity, etc., then it is usually 
more practical to use the Kadanoff-Baym equation and its ensuing Boltzmann- 
type approximations to be studied in the next sections [16]. Then one has to 
solve an (integro-) differential equation. It is easier to make approximations 
(such as near-homogeneous or slowly-varying perturbations, not too far from 
equilibrium) on the Kadanoff-Baym equation than on the Keldysh equation. 
Furthermore, a connection to conventional transport formalisms such as linear 
response theory (Kubo formula) can be made: transport coefficients calculated 
from the linearized quantum Boltzmann equation (which can be derived from 
the Kadanoff-Baym equation) in the steady-state limit are identical to those 
obtained in linear response theory. 

2.2 Kadanoff-Baym Equation 

With the help of these relations, we are now able to derive the Kadanoff-Baym 
equation from the Keldysh equations. If we act with Gq 1 on the left of the 
first Keldysh equation and on the right of the second Keldysh equation and 
use relations (2.15)-(2.18), we obtain 

{Gq 1 - U)G< = U R G < + E<G A (2.21) 
G < (Gq 1 — U) — G R S < + G < S A (2.22) 

We now subtract these equations from one another: 

[Gq 1 - U, G<] = E R G< + U<G A - G R E< - G<E A (2.23) 

We now introduce the nonequilibrium spectral function A = i(G R — G A ) 
and the scattering rate, or linewidth, or imaginary part of the self-energy 
r = i(E R — S A ). We also define symmetric combinations or real parts as 
ReG = \{G R + G A ) and ReU = \{S R + U A ). In terms of these quantities, 
Eq. (2.23) can be written as 

[Gq 1 - U, G<] - [Re 27, G<] - [£<, ReG] = -\ ({r, G<} - {A, £<}) 

However, as seen before, the Keldysh components G R,A,<,> and S R,A,<,> 
are not independent. With the relations — iA — G R — G A — G > — G K and 
— iF — S R — S A — — S K we can rewrite the equation as 

[Gq 1 - U -ReU,G<} - [S<,ReG] = \ ({E>, G<} - {G>, (2.24) 

which is the Kadanoff-Baym equation [8, 13]. As we will see, the term on 
the left-hand side involving [Gq 1 — U, G K ] is a driving term, while the terms 
involving ReU and ReG describe renormalization effects. From a transport 
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point of view, they lead to renormalized transport coefficients and are ne- 
glected in the classical Boltzmann limit. The right-hand side is a collision 
term as will be seen. A similar Kadanoff-Baym equation can be derived for 
G > . By subtracting the two equations, one obtains an equation for the spectral 
function, 

[Gq 1 - U-ReU,A] - [r,ReG] =0 (2.25) 

which is used as a consistency check when one looks for approximate solutions 
to the Kadanoff-Baym equation. 

The Kadanoff-Baym equation has to be supplemented by an equation for 
the retarded and advanced Green's functions. It is obtained in the same way, 
applying relations (2.15)-(2.18) to the nonequilibrium Dyson equations (2.3)- 
(2.4): 

(Go 1 - U)G R = 1 + E R G R (2.26) 
G r {Gq 1 -U) = l + G R U R (2.27) 

Add these two equations and dividing by 2 gives 

\{Gq 1 - U, G R } = 1 + \{Z R , G R } (2.28) 

In practice, because self-energies are typically functionals of the Green's func- 
tion U — £[G], the Kadanoff-Baym equation and the equations of motion 
(2.26)-(2.27) for G R have to be solved self-consistently. Needless to say, this is 
a rather difficult task and one often has to resort to numerical techniques. Ap- 
proximations have to be made for the self-energy functions. However, these 
approximations cannot be made blindly by choosing an arbitrary subset of 
diagrams since it is then possible to violate conservation laws as shown by 
Baym and Kadanoff. Self-energies should be derived from the Luttinger-Ward 
[20] functional 0, 

m = ^ (2-29) 

The Luttinger-Ward functional $[G] is defined diagrammatically as the sum 
of all skeleton connected vacuum diagrams with free propagators Go replaced 
by exact propagators G. In practice, one keeps only a few diagrams in the 
diagrammatic expansion of and then obtains the corresponding diagram- 
matic representation of the self-energy S[G] by removing a fermion line G 
(which corresponds to functional differentiation in Eq. (2.29)). These so-called 
^-derivable approximations preserve conservation laws and are termed con- 
serving approximations [21, 22, 8]. For interacting systems involving coupled 
degrees of freedom such as electrons and phonons, the electron self-energy U 
and phonon polarization 77 should be derived from the same electron-phonon 
vertex r. This ensures that the same level of approximation is maintained 
for both the electron and phonon subsystems [23]. In general, gauge invari- 
ance implies the Ward identities which relate the self-energies and the vertex 
so that a given ^-derivable approximation for the self-energy U[G] defines a 
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fully consistent level of approximation for Green's functions, self-energies, and 
the vertex. 



2.2.1 Wigner Representation and Gradient Expansion 

We now move on to the gradient expansion of the Kadanoff-Baym equation 
[8, 16, 24]. The idea is to separate slow macroscopic variations from fast 
microscopic variations, and then perform a gradient expansion on the slow 
variables. In order to do that, we introduce mixed coordinates, the so-called 
Wigner representation: 

r = xi-x 2 , R=i(xi+x 2 ) (2.30) 
t = t!-t 2 , T=|(ti+t 2 ) (2.31) 

so that a function C(xi, t\\ x 2 , t 2 ) becomes a function C(r,t,R, T). The fast 
'relative' variables (r,t) will be Fourier transformed to (k, Q) while the slow 
'center-of-mass' variables 3 (R,T) will serve for the gradient expansion. The 
Fourier transform is defined as 

C(k,X2,R,T) = JdtJ dre^- k - r) C(r,t,R,T) (2.32) 

= J dt J (ire i( ^- k - r) C(R+|r,r+^;R-|r,r-|t) 

We first work out the Wigner transform of the driving commutator term 
[Gq 1 - U,G<}. We have 

[Go 1 - U, G<] X1 , X2 = (i J- - fT (xi, -*Vi) - U( Xl )j G < (x 1 ,x 2 ) 

-G < (x 1 ,x 2 ) ~ ^o(x 2 ,«V 2 ) - U(x 2 )J 

and for convenience we use a four-vector notation x± 2 = (ti, 2 ,xi ;2 ). 

We want to study the dynamics of a system driven out of equilibrium 
by constant and uniform external electric E and magnetic B fields 4 . We first 
concentrate on the electric field and the magnetic field will be added later. The 



3 For example, consider a free propagator g(t±,t2) = — ie — ' iro< -* 1— * 2 - ) where zu is 
some characteristic frequency. Considered as a function of t = |(ti — £2) and 
T = |(£i + £2), this function has a fast variation with frequency 2zu in the t 
variable and is constant (infinitely slow variation) in the T variable. If we add a 
slow disturbance vj — * w + A(t), there will start to be a weak dependence in T, 
but slower than the dependence in t. 

4 All the following derivations can be done for a general one-particle potential U(x) 
but would lead to extremely cumbersome expressions. This doesn't have much 
interest anyway since the physically relevant perturbations are E and B fields. 
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one-particle potential U{x) contains the electric field in the scalar potential 
gauge: 

U(x) = -eE-x (2.33) 

We have 

[Gq 1 -U,G<] X1 , X2 = I i— +z — J G < (x 1 ,x 2 ) - H {K 1 ,-iV 1 )G < {x 1 ,x 2 ) 

+G < (xi,x 2 )-£/o(x2,^V2) + eE ' ( x i - x 2 )G < (xi, x 2 ) 

The change of variables Eqs. (2.30) and (2.31) imply the following relations 
for the derivatives, 

d d d d 1 / d d 

and similar relations for the spatial derivatives. In Wigner coordinates we thus 
have 

[Go 1 - U, G < ] r>t> R, T = ^ G< ( r > *> R > T ) ~ H o (R + |r, -*(§V R + V r )) G< 
+G < # (R - §r, i(|V R - V r )) + eE • rG< 

We now perform the Fourier transformation of the center-of-mass variables 
Eq. (2.32), which corresponds to the following substitutions, 

d 

i— -> Q (2.34) 
dt v 7 

-zV r -> k (2.35) 
t^-iA (2.36) 
r -> zV k (2.37) 

We obtain 

[Go 1 - 17, G<} Kn ,n,T = i ^G < (k, Q, R, T) -H Q (R+ |V k , k - |V R ) G< 

+G<i/o (R - f V k , k + | V R ) + zeE • V k G< 

We now perform a change of variables introduced by Mahan and Hansch, 

w = i? + eE-R (2.38) 

which eliminates an unphysical oc E • R term in the equation for G R to be 
derived later. This change of variables implies the following substitutions, 

Q -> u - eE • R (2.39) 
V R ^V R + eE-|- (2.40) 
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The derivative d/doa is important and differentiates the QBE from the classical 
Boltzmann equation. We thus obtain the following final exact form for the 
driving commutator in the Wigner representation, 

[Go 1 ~ U, C < ] k ,.,R, T = - H (R + §V k , k - | (V R + eE^-) ) G< 

+G<fr (r - | V k , k + | (|vr + eE^- J J + zeE • V k G< 

(2.41) 

where G < = G < (k, u, R, T). We now see the importance of a proper ordering 
of the Hamiltonian since [R, Vr] ^ and [k, V k ] ^ 0. The driving anti- 
commutator {Gq 1 — U, G R } appearing in Eq. (2.28) for the retarded Green's 
function can be derived in a similar way. We simply quote the result: 

IIGo 1 " U, G r } Kuj ^t = ojG r - \Ho (r + §V k , k - | (v R + eE^ ) G R 



\G R H I R - §V k , k + | ( V R + eE^ ] ] (2.42) 



d_ 



where as advertised the unphysical term E • R is removed by the Mahan- 
Hansch transformation Eq. (2.38). 

Having derived the driving terms in the Wigner representation, we must 
now see how to transform products of the form AB, 

(AB) Xl , ti;X2 , t2 = J dt' J dx' A(x 1 ,t 1 ;x',t')B(x',t';x2,t 2 ) (2.43) 

From the definition of the Wigner coordinates Eqs. (2.30-2.31), we have 
A(xi,ti;x',t') = - x',t! - if, |(xi + x'), \{ti+t')) 

= ^(x 1 -x^ 1 -^R+ 3C ^,T+^) (2.44) 

and similarly 

B(x', t'- x 2 , t 2 ) = 8(x' - x 2 , t'-t 2 , §(x' + x 2 ), §(*' + t 2 )) 

= s(x'-x 2 ,t'-t 2 ,R+^^,T+^) (2.45) 

Since the derivative is the generator of translations, we have the following 
compact form for the usual Taylor expansion, 

/(R + a, T + s) = e a - VR e sdT f{TL 7 T) (2.46) 

Using Eqs. (2.44,2.45,2.46), we can write Eq. (2.43) as 
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(AB) Xl , ti;X2 , t2 = Jdt' J dxM( Xl -x',*! -t',R,r)e-^- x > v -/ 2 e-^-^/ 2 
xe (x'-x 2 ).v^/2 e (t'-t 2 )a#/2^ (x / _ t f _ ^ R r) (2>47) 

where the superscripts A and B on the derivative operators Vr,<9t indicate 
that they act on the function A and B, respectively. The integral now has the 
form of a convolution product, 

J dt' J dx! I(xi - x', ti - t', R, T)£(x' - x 2 , £' - £ 2 , R, T) = A* B, 

where 

I(r, R, T) = A(r, t, R, T) e - (r - v *+^/ 2 (2.48) 
£(r, t, R, T) = e ( r - v £+^)/ 2 £(r, t, R, T) (2.49) 

We now take the Fourier transform, as defined in Eq. (2.32), of the convolution 
product, so that we have 

(AB) k ,^ R)T = T{A * B} = A{k, n, R, T)£(k, f2, R, T) (2.50) 

since the Fourier transform of a convolution product is the ordinary product 
of the Fourier transforms of the convoluted functions. We thus have 

A(k,n,K,T) = J dt J dre 2(m - k - r U(r,t,R,r)e-( r - v ^+^ ) / 2 (2.51) 

B(k, n, R, T) = J dt J dr e Wt-*- r ) e ( r -v& +td 2V 2 B(r, t, R, T) (2.52) 

Since a phase factor in a Fourier transform generates a translation, 

F{e a - r e st f(r, t)} = /(k + ia, Q - is) = e ^V k£r ^ f( kj ^), 
we obtain 

I(k, a, r, r) = A(k, a, r, T) e -* v *' v * /V^ /2 (2.53) 

B{k,n,K,T) = B{k,f2,-R,T)e iV ^E/2 e -^d? 2 /2 ( 2 .54) 

Inserting these results in Eq. (2.50), we have 

(AB) k ,rt,R, T = A(k, 12, R, T)£ AB (k, 12, R, T)B(k, ft, R, T) (2.55) 

where we have now dropped the curly letter notation A and B for simplicity, 
and we define the gradient operator Q AB as 

Q AB {k, H, R, T) = e - i ( 9 T^-^9^-v^-v^+v^-v|)/2 

As a mathematical aside, we note in passing that Eq. (2.55) corresponds to 
the Moyal or star product of two functions which is used in noncommutative 
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geometry and deformation quantization, and more generally in mathematics 
in the construction of deformed algebras. 

So far, everything has been exact. In the so-called gradient approximation, 
we expand the exponential to first order in the gradients. We have 



e - (k)12)R)T) ^_«(|_|__|_|__ v .. VkB + v .. vi 



so the basic expression of the gradient approximation is 



(Afl) k>fl>R> T A(k, n 7 R, T)B(k, R, T) 
i fdAdB 8A8B 



2 \dT dQ dQ dT 



V R A • V k E + V k A • V R 5 



from Eq. (2.55). We can now write down the Kadanoff-Baym commutators 
and anticommutators in the Wigner representation under the gradient ap- 
proximation as follows, 



[A,B]^,t = [A,B]+ 1 -(^ 



+ -eE 



{A,B} Ku) ,h,t = {A,B} + 



+ -eE 



\dA dB 

aJ' dT 

OA 



OA dB 

d ^ B \-\^ A ' 9 £ 





'OA 


dB' 




'OA 


dB' 


( 


doo ' 


df_ 






duo 



OA 



,V k E 



dB 



-{V k AV R B} + {V R i,V k B} 

(2.56) 

[V k A,V R 5] + [V R A,V k B] 

(2.57) 



where on the right-hand side, these are ordinary matrix commutators and 
anticommutators, and we have performed the Mahan- Hansen transformation 
Eq. (2.38). 



2.3 Quantum Boltzmann Equation 

We first discuss the QBE [24] in zero magnetic field, then add the magnetic 
field. 



2.3.1 QBE with Electric Field 

The QBE is now obtained merely as the Kadanoff-Baym equation (2.24) in 
the Wigner representation: 

[G — U, G < ] k)U ; i R j T = [Re E, G < ] k + [U<,ReG\ k ,u>,R.,T 

+ 1 ({^' > ' ( ^ < }k,u;,R,T - {G > , X' < } k ,^,R,T)(2.58) 
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where the left-hand side is given by Eq. (2.41) and the terms on the right-hand 
side are given by Eqs. (2.56) and (2.57). So far we have only rephrased the 
Kadanoff-Baym equation in Wigner coordinates and made the gradient ap- 
proximation. Equation (2.58) is still valid to arbitrary order in the electric field 
and for both time and space-dependent (but slowly varying) perturbations. 
Note that this very general equation has in principle to be solved together 
with the equation for the retarded Green's function Eq. (2.28), 

H^o 1 ~~ U, G r }\s_^^i,t — 1 + \{^ R i G r }\l^,yi,t (2.59) 

where the left-hand side is given in Eq. (2.42) and the anticommutator on the 
right-hand side is given by Eq. (2.57). Indeed, the nonequilibrium G R enters 
the renormalization terms KeG and KeU (since U R — U R [G R ]) in Eq. (2.58). 

It is obvious that this primary form of the QBE is extremely complicated 
and that not much progress can be made unless further approximations are 
introduced, otherwise one has to resort to numerical techniques. In addition, 
we have kept a fully general multiband Hamiltonian Ho so far, but one can 
also restrict the analysis to a given class of Hamiltonians. 

2.3.2 QBE with Electric and Magnetic Field 

We now add the effect of a constant uniform magnetic field B. Whereas the 
electric field E was introduced in Eq. (2.33) through the scalar potential, we 
now introduce the magnetic field through a vector potential, 

A(x) = -|x x B 

We also perform the Peierls substitution in the Hamiltonian, 

H (-iV) -> H (-iV - eA) = H (-zV + §ex x B) (2.60) 

We now proceed along similar steps as before. In Wigner coordinates, we have 

H (x u -iVi - eAi) -> H (R + |r, -i( JVr + V r ) + §e(R + |r) x B) 
H {x 2 ,iV 2 - eA 2 ) -> # (R - |r, i(|V R - V P ) + |e(R - §r) x B) 
After Fourier transformation, we obtain terms like 

H (R ± | V k , k + |eR x B t |(V R + \eB x V k )) 

As before, we perform the Mahan-Hansch transformation Eq. (2.38) which 
changes the derivative with respect to R according to Eq. (2.40), and obtain 

H (r ± f V k , k + §eR x B =f f (^V R + eE-|- + |eB x V k 

The Mahan-Hansch transformation got rid of the oc E • R term. We now 
perform a similar extra transformation to get rid of the unphysical R x B 
term by introducing the kinematical momentum, 



) 
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p = k+ |eR x B = k- eA(R) (2.61) 
which implies Vk — > V p but modifies the R derivative once more, 

Vr^Vr+ |eB x V p 
so that the final expression reads 



fin [ R 1: ^ V p , p =F ^ j V R + eE-|- + eB x V p 



Hence the generalization of Eq. (2.41) to include the magnetic field B is 

[Go 1 - U, G<] p ^r,t = i 9 ^- - H (k + | V p , p - | (V R + eE^- + eB x V p ) ) 



+G<H ( R - | V p , p + \ ( V R + eE^- + eBxV, 



+zeE • V P G < (2.62) 
and similarly the generalization of Eq. (2.42) is 

\{G^ - U, G R } P ^ R , T = ujG r - \H Q (r + | V p , p - § (V R + eE^- + eB x V p ) ) G K 

- \G R H (it - | V p , p + f fv R + eE^- + eB x V p ) J (2.63) 



All Green's functions and self-energies are now functions of (p, uj, R, T). The 
last thing to figure out is how Eqs. (2.56) and (2.57) are modified in the 
presence of the magnetic field. It is straightforward to show that the only 
changes are the substitution k — > p and a new term linear in the magnetic 
field, 

([A,fi]±)p, w>R>T = ([A,fi]±) k>w>R , T + |eB • [V P A x V p £]± (2.64) 

where we use the notation 

[V P A x V pJ B]± = V P A x V pJ B ± V P B x V P A (2.65) 

Then the QBE in the presence of both electric E and magnetic B fields is 
formally identical to Eq. (2.58), 

[G Q 1 — U, G < ] p a;j R 5 T = [Re S, G < ] PiU ; j r 5 t + [^ ,< , Re G] PjLJ ,n,T 

+ | ({i7>,G < }p >w> R,r - {G>,r<}p iC ,,R T j2.66) 

but with the definitions in Eqs. (2.62) and (2.64). Similarly, the equation for 
the retarded Green's function Eq. (2.59) becomes 
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|{^0 —U 7 G }p jU ),R.,T — 1 + \{£ ,G } p ,^,r,t 



(2.67) 



together with Eqs. (2.63) and (2.64). 

Equation (2.66) is the most general form of the quantum Boltzmann equa- 
tion and makes no assumptions other than that of slowly varying disturbances. 
It is obvious that it is a rather complicated equation. Approximations should 
be introduced to bring it to a more tractable form. The first obvious approx- 
imation is to restrict the analysis to linear transport, that is, keep only terms 
to linear order in the electric field, while keeping terms to all orders in the 
magnetic field. Afterwards, approximations can branch off in different direc- 
tions and depend on the problem at hand. In the next chapter, we discuss 
several useful approximate forms of the QBE. 



2.3.3 One-Band Spinless Electrons 

Most analytical work using the QBE is done for one-band spinless electrons. 
The single-particle Hamiltonian Hq is of course 



tfoHV) = 



HV) ; 
2m 



The driving terms Eqs. (2.62) and (2.63) become 

[G Q 1 — U, G < ]p, w ,R,T = 
|{^0 1 ~~ G' R }p,a;,R,T — 



d_ 
df 



d_ 



+ v p • I V R + eE— ) + e(E + v p x B) • V p 



1 / „ d 



eB x V, 



G 



where we have the usual definitions of the single-particle energies and velocities 

.2 _ 



£n = 



P 

2m 



v p = V p e p = 



rn 



We see that the driving term [Gq 1 — U, G K ] has indeed the form of Eq. (2.2) 
derived from simple considerations. We now consider the Kadanoff-Baym com- 
mutators Eq. (2.64). Consider first the anticommutator {A, 5} p ,u;,r,t- All the 
commutators in Eq. (2.57) vanish since the quantities are all scalars. Further- 
more, the term proportional to B in Eq. (2.64) vanishes because of the anti- 
symmetry of the cross product. Hence we have simply {A, -B} p , w ,r,t = 2AB, 
and Eq. (2.67) becomes 



w - e D + 



8m 



V R + eE^- + eB x V p 



-S 1 



G R = l 



(2.68) 



The Kadanoff-Baym commutator in Eq. (2.64) becomes 
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r„ ^ fdAOB OA OB „ „ „ „ „ „ „ \ 

+zeE • f ^jV pJ B - V P A— J + zeB • (V P A x V P B) 

so that the QBE Eq. (2.66) becomes 

dReU\ d dReU d „ n n , / „ 9 

+ (v p + VpRei:)- V R + eE- 



duj J dT dT Ouj v p p y V du; 
+e^l-^^jE + (v p + VpRei;)xBj • V p - V R Re T • V p G< 

f p)y< f)~RpG \ 
-z'eE • f V p Re G — V p r< j - ieB ■ (V p r< x V p Re G) 

where we use the four- vector notation p M = (oo, p) and X M = (T, R) with the 

Minkowski metric r]^ v — (H ). 

Equation (2.69) is valid to arbitrary order in E and B fields. We see that 
the self-energy term [Re U, G K ] in the right-hand side of Eq. (2.66) has renor- 
malized the velocity in the driving term, 

v p -> v p + V p Re E = V p (e p + Re S), 

and introduced the factor (l — d ^ S ) reminiscent of the wavefunction renor- 
malization factor Z _1 in Fermi liquid theory. The other term [X' < ,ReG] 
gives some additional terms proportional to the E and B fields. The term 
E > G < — G > U < on the left-hand side is a collision term. 



3 



Applications 



We now consider a few applications of the nonequilibrium formalism. 

3.1 Nonequilibrium Transport through a Quantum Dot 

In this section we study nonequilibrium transport through a quantum dot 
connected to two external metallic leads (two-probe system). Let us first derive 
a general expression for the current through the dot in terms of Keldysh 
nonequilibrium Green's functions. 

The Hamiltonian of the whole two-probe system is 

ka£L,R ka£L,R,n 

(3.1) 

where c^ a (ck«) is a fermionic creation (annihilation) operator for a single- 
particle momentum state k in channel a in the left or right metallic lead, and 
d^jdn are creation/annihilation operators for states in the quantum dot. He 
is the Hamiltonian of the quantum dot. For a noninter acting dot, Hq would 
be 

H C = ^2 h"mndtnd n (3.2) 
mn 

For a dot with on-site repulsive interactions and a single site with spin-split 
levels e a , we use the Anderson model, 

H c = ^2e a d^ a d a + Un^ni (3.3) 

a 

with n a — d\d a is the number operator on the dot. 

In any case, let us first derive a general expression for the current regardless 
of the presence of interactions in the dot. We however do require that the 
interactions, if any, be limited to the electrons inside the dot: the leads should 
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be noninter acting and there should be no interactions between electrons in 
the leads and electrons in the dot. 

The full Hamiltonian is partitioned as explained in Chapter 1, i.e. Ti — 
H + H' where H is the equilibrium Hamiltonian and H' is the nonequilibrium 
perturbation. We also have to specify the initial noninteracting density matrix 
p{— oo). In the present problem, we have two possible ways of doing this 
partitioning: 

• Choose H to be the Hamiltonian of the isolated leads and dot, and H' to 
be the coupling between the dots. The initial density matrix in the remote 
past p{—oo) is the product of the density matrices of the isolated systems 
in equilibrium: 

H = ^2 e kaCk a c ka + H c {{dl}, {d n }} 

ka<EL,R 
kaEL,J2,n 

P(-OO) = 1 (e~^-^N L ) % p (0) ^ e -(3 { H R -» R N R )^ (3 4) 

where Hl,r are the Hamiltonians of the isolated leads in equilibrium each 
with chemical potential Pl,r, and the initial density matrix of the dot 1 
p^) involves only the noninteracting part of the Hamiltonian of the dot 
H c - 

• Alternatively, we can choose H to be the Hamiltonian of the connected 
system in equilibrium at a single chemical potential p. The nonequilibrium 
perturbation H' is then a one-body potential term which raises the single- 
particle energies by in the leads. This represents the shift 2 in chemical 
potential p — > Pl,r = P + ^l,r- We then have 



1 The obvious choice is p£ = e~ l3<yIi c ~^c N c) w h ere jj^> i s the noninteracting 

< 

part of He but pc is unspecified. However this initial condition is included in Gq 
which drops out of the Keldysh equation, see Section 2.1, so that the steady-state 
at any finite time t > — oo after the Keldysh adiabatic evolution is independent of 
pc- In particular, the steady-state current is independent of pc- This is a example 
of washing out of the initial conditions in the Keldysh formalism (even if here 
we are not even considering initial correlations) . The steady-state nonequilibrium 
population of the dot will be given by the nonequilibrium Green's function of the 
dot G<. 

2 To get the correct nonequilibrium self- consistent potential profile throughout the 
device including the dot, one would have to include Coulomb interactions in the 
Hamiltonian, at least at the Hartree (mean-field) level. Otherwise one can neglect 
interactions and add a term in the nonequilibrium perturbation to mimick the 
self-consistent potential profile, 5H' = Ad) d where A = Al + Ar ; assuming a 
symmetric device. 
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H = e kaC^ a C kce + Hc[{dl}, {d n }] + ^ \t^a,nC^ a d n + ^ka,n^n C ka 

kaGL,i? 

p(-oo) = \ e -m-^) 

We will choose the first method which is well-suited to steady-state problems, 
while the second method is easier for time-dependent transport. 

3.1.1 Expression for the Current 

The current flowing through lead a is defined as 

J a (t) = -e(N a ), (3.5) 

where N a — X]k c ka Ck « * s ^ e number operator for lead a and e > is the 
electron charge. From the Heisenberg equation of motion ifiN a — [N a ,H] and 
Eq. (3.1), we obtain 

J a (t) = ^Re^t ka , n G< ka (t,t), (3.6) 

k,n 

where we have the mixed lesser Green's function G^ ka (t, t') — i{(^_ a {t')d n {t)) . 
Using the Keldysh technique, we will obtain an expression for the mixed 
contour-ordered Green's function G nika (r, t') — — «(r c {d n (r)c kQ ,(r / )}) and 
perform analytic continuation to real time to obtain the lesser function 

G< ka (M')- ' 

For pedagogical reasons, let us first derive the expression from the usual 
perturbation expansion of the S'-matrix and the subsequent application of 
Wick's theorem, and then obtain it from the path integral method which is 
more transparent. 

3.1.2 Perturbation Expansion for the Mixed Green's Function 

As explained, we define the unperturbed Hamiltonian as the first two terms 
of Eq. (3.1), 

H = E £k " C ka C ka + H C [{dl}, {d n }], 
ka 

and the perturbation as the tunneling term, 

H — ^ ^ (^t'ka,ncl ( _ a d n + £ kQ!)n ^nCka^ • 
ka,n 

The perturbation expansion becomes 
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G n ,Ur,T , ) = J2- Z:! Tr- f dn---f ^(T c {d n (r)^(r 1 )---^(r z )cL(r / )}), 
1=0 ' ^ c ^ c 

(3.7) 

where we have 

i 

(T c {d n (r)H'(r 1 )---H'(r l )cl a (r')}) = [I E ( T ° (4<«<fa)*i*«*,»A 4 (T< 

i—l kiai ,rii 
+4 i (Tt)*k i a i ,nAi«i( T i)) 4a ( r ')}) • 

Wick's theorem can now be applied to correlators of c, & fields since if is 
quadratic in these fields, whence the importance of the requirement that the 
leads be noninter acting. All the contractions involving the first term in H'fa) 

n, n, 

vanish, since d n c\^. a . — and d ni c\^ a — for the unconnected system, and the 

anomalous contractions (inrfn, and c^ ^.c^ vanish because particle number is 
conserved. By Wick's theorem we thus have 

(T c {4(r)^(r 1 )... J ff'(^)cLM}> 
z 

= 11 E ^^(^{^(^^(rOc^a^rOaLM}) 

+ E *k,_ 1 a l _ 1> n,_ 1 (T c {4(r)^(n) • • • H'in^H'in)^ (r^)}) 

kj_ia;_i,nj_i 

x<T c {4 I _ iai _ 1 (7 l _ 1 )cl Ea (r / )}> + ... 

kiai ,ni 

where all the I terms are seen to yield the same contribution to the pertur- 
bation expansion Eq. (3.7) if the dummy integration variables T\, . . . , t\ are 
relabeled. Choosing the labeling of the last term, we obtain 

C n ,ka(r, t') = E £ dn(-i) Jj £ <Zt 2 • • • £ drz 



c 



After relabeling the dummy integration variables once more, we see that the 
perturbation terms H' ■ ■ ■ H' give rise to the contour ^-matrix Eq. (1.26), so 
that we finally have 

(n,r), (3.8) 
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where 

G nm (r,r') = -i(T c {d n (r)dl(r')}} = -i(T c {S c (-oo, -oo)d n (r)dl(r')}), 

(3.9) 

is the contour-ordered Green's function of the scattering region, and we have 
-i{T c {cK a {T)c^, al {T')}) = <W<W0ka(T,r'), 

because of translational invariance in the leads, where g^ a (T, t') = — ^(2c{cka(T)c^. Q ,(r / )}) 
is the contour-ordered Green's function of the (isolated) leads. 

3.1.3 Path Integral Derivation of the Mixed Green's Function 

The path integral can also be used in the Keldysh technique. The only dif- 
ference is that the action is obtained by integrating the Lagrangian over the 
Schwinger-Keldysh contour C instead of over the real axis. The Lagrangian is 

L(c, c,d, d) = y^c ka (id T -e ka )c^ a +y^d n id T d n -H c [{d n }, {d n }]~ (cka^ka 

,ndn~\~d n t-^ an C]^ cx ^ 

ka n ka,n 

(3.10) 

To calculate the propagator G n ^ a (T,T'), it is appropriate to define the 
Keldysh generating functional Z[fj, J], 



Z% J]=TrpT c 



where f\ n and Jk a are Grassmann sources, and the normalization is such that 
Z[0, 0] = Tr p = Z. The generating functional is useful because it can be used 
to generate the mixed correlation function, 



iG n) k a (r,T') 



1 S'Z%J] 



(3.11) 

77=0, J=0 



The generating functional has a path integral representation, 

Z% J} = J ^[c,c,d,d]^( S[ "' c ' J ' d]+ i dT(??d+ " J) ) 

where we use the shorthand notation fjd = ^ n f\ n d n and cJ = Xlka ^ka^ka- 

The d and the c fields are entangled in the Lagrangian Eq. (3.10). To 
disentangle them, we perform a shift of variables, 

4a( T ) = Cka(r) ~J2f dn dm(Ti)££ a m 0ka(Tl,T) 
m J C 

C ka( r ) = Cka(r) dT ± 9ka(T, T\ )tka,m^n fa ) 
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where <7ka(r, T\) is the contour-ordered Green's function of the isolated leads 
defined earlier. It satisfies (see Eqs. (2.9) and (2.10) of the notes) 

(id T - e ka )# ka (r, Ti) = 5{t, Ti) 
#ka(r,Ti)(-z<9 Tl - e ka ) = 5(t,ti) 

where in the second equation the derivative acts to the left. With this shift of 
variables, the Lagrangian becomes 

L(c,c',d, d) = V"c k a(^r-eka)c ka +L ({(in}, {d n })- V" f dTid m (r)tl am g^ a (T, Ti)£ ka , n d n (ri) 
, i Jc 

where L ({d n }, {d n }) = Y, n d n id T d n - H c [{d n }, {d n }] is the Lagrangian of 
the isolated dot. 

EXERCISE. Show this (hint: use integration by parts). 
As a result, the action becomes 

S[c', c', d, d] = S^ads [c', d] + S QD [d, d] 

where 



Sleads [cV] = yZf drc '^ a { id r ~ £ka)c ka 
ka J C 



is the action of the isolated leads, and 



S QB [d 7 d} = f dr L ({d n } 7 {d n })- V * dr f dr x d m (r)£ kam # ka (T, ri)t ka>n d n (ri] 
Jc , — Jc Jc 

ka,mn 



is an effective action for the quantum dot, which contains the effect of the 
external leads (bath). The source term f)d is unaffected by the change of 
variables, but the source term cJ becomes 

CJ = ^C< ka {T)J^ 0i {T) = I 4a ( T ) + f rfr l^m(Tl)t ka)m ^ k a(Tl,r) I J ka (r) 

ka ka V m J ° J 

The generating functional thus becomes 

Z% J}= J V[c\c\d,d]e^^'^e tS ^ d ^e 

where the integration measure T>[c, c] —T>[c', d] is invariant under the change 
of variables since it is only a shift. 

We now perform the functional differentiation Eq. (3.11). This generates 
two terms. The first term is 
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\ J ^[c^ C ^J,rf]e^e^^4r)4JrO = (^Jd n (r)4 a (rO}) 

which is zero because it corresponds to the unconnected system (the action 
in the path integral is that of isolated leads). The second term is 

m 



where 



G„ ro (r,n) = <T c {d n (T)^(n)}) = iy 2?[fi',c',J,d|e <Sl — e^^^^Cn) 
£>[c, c, d, d]e lS d n (T)d m (Ti) 
is the exact Green's function of the connected dot. Hence we obtain 
(r, t) = \ & dri G nm (r, Ti)t£ a (ri,r') 

m ^ 

which is just Eq. (3.8). 

3.1.4 General Expression for the Current 

In Eq. (3.6) for the current, we need the lesser mixed Green's function. We 
therefore apply the Langreth analytic continuation theorem Eq. (1.64) to Eq. 
(3.8) to get 

G n ^ a (t,t') = / dt 1 (G^ m (t,ti)tl ajm g^ a (t 1 ,t')+G^ m (t,ti)tl aim g^ a (ti,t')) 

m 

We now assume steady state so that all the Green's functions depend only on 
the time difference, i.e. G(t,t') — G(t — t'). Then we get 

Gn,ka(M)=E / ^(Cl^^o,^^) + G< ra ( W )tL,™rfa(^)) (3-12) 
m 

The Green's functions of the isolated leads are 

0kaM = 2TTif a {e^ a )5(uj - e ka ) (3.13) 

9tH = ; * ; (3-14) 
uj - id - e ka 

where f a (e) — (e^( e_Ata ) + 1) _1 and the chemical potentials are those of the 
initial density matrix Eq. (3.4). We now substitute Eq. (3.12) in Eq. (3.6), 
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J a = yRe^ y'^tL,mW,n(PkaM G ?mM+^MG ! nmM) ( 3 - 15 ) 
k,mn 

It is not hard to show from the definition of the Green's functions that the 
following relations hold, 

C < H t = -G<{uj) 
G R (uj)i = G a {uj) 

considering the Green's functions as matrices G — G nm . Using these relations 
and Eqs. (3.13) and (3.14), we can show that Eq. (3.15) becomes [25] 

Ja = j J ^Tvr a {u){Uu)[G R {u)-G A {u)]+G<{u)} (3.16) 

where we define a linewidth function 

and we have introduced the density of states p a {w) = ^ k 5(lv — eka) to convert 
the sum over k in Eq. (3.15) to an integral, 

J2 F ( 6 ^) = [ de Pa (e)F(e) 
k ^ 

Equation (3.16) is an exact expression for the steady-state current through 
lead a of an arbitrary multiprobe system with interactions inside the dot (but 
not in the leads). However, the nonequilibrium Green's functions have to be 
calculated in the presence of interactions, and in the presence of the leads held 
at different chemical potentials. 

For a two- probe system with proportionate couplings to the leads, -Tl(^) = 
A.Tft(u;), one can arrive at a simpler expression. We first use the fact that in 
steady state, J — Jl — —Jr, so that we can write J = xJl — (1 — x) Jr for 
an arbitrary x. The current then reads 

J=^J ^Trr R [(\x-(l-x))G< + (\xf L -(l-x)f R )(G R -G A )] 

We then fix the arbitrary parameter x — so that the first term vanishes 
and the current does not depend on G< anymore. We then obtain 

' " ~n I >^ ~ '»«"» * (tStS) ) ^ - G »i ^ 

where the ratio is well-defined since the matrices F^ and Fr were assumed to 
be proportional. 
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3.1.5 Noninteracting Quantum Dot 

We derive an alternate expression for the current in the case that there are 
no interactions in the dot. In this case, the retarded Green's function satisfies 
the Dyson equation 

G R = G R + G R S R G R 7 (3.18) 

where 

^nM = E^«,^kaMW (3.19) 
ka 

is the noninteracting tunneling self-energy, with r — ^pTp — i(U R — U A ). 
G R is the equilibrium Green's function of the unconnected dot, 

G r (uj) = ((v + iS-h)- 1 . 

In other words, in the noninteracting case, the self-energy contains only the 
effect of the leads. From these equations it is easy to derive the identity, 

G R — G A — -iG R FG A 

The lesser Green's function follows from the Keldysh equation, 

G< = G R U<G A 

where the lesser self-energy is 

^^E^ ( 3 - 2 °) 

Consider a system with a single level, such that all quantities G, r, U are 
scalars. If we write G< in pseudoequilibrium form 

G< = if A 

with A — i(G R — G A ) — G R FG A the spectral function, the pseudodistribution 
/is 

j _ fh^L + [rTr 

1 ~ r L + r R 

which is clearly not of equilibrium form if Jl ^ Jr. 
Using the previous results, one can show that 

J ot = V^ I duj[f a {uj) - ff3{uj)]T a p{uj) 

P 

where the transmission coefficient T a/ 3 is 

T a/3 (u;) = Trr a (u)G R (cj)r fi (u)G A (u) 

which is the usual Landauer formula. Note that Landauer obtained this result 
heuristically while here it has been derived rigorously from the nonequilibrium 
formalism. Also, the Green's functions G R and G A have to be calculated out 
of equilibrium. 
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3.1.6 Interacting Quantum Dot: Anderson Model and Coulomb 
Blockade 

In this section we consider on-site repulsive interactions in the dot [26, 27]. We 
assume that the couplings are proportionate F L oc so we use Eq. (3.17). 
The Hamiltonian of the quantum dot is the Anderson model Eq. (3.3), 

a 

Assuming that the couplings are diagonal in spin space, the current is 

J = lE/ MU(uj) - /«M]r w M (-iimG^)) 

where 

G* a (t) = -ie(t)({d a (t),dUo)}) 

is the nonequilibrium retarded Green's function of the dot (in the presence of 
the leads). This Green's function also satisfies a Dyson equation Eq. (3.18), 
but the self-energy now contains both contributions from the leads and from 
interactions in the dot. Both contributions are entangled in the exact self- 
energy which is not known. Obviously because of interactions the problem 
cannot be solved exactly. We will solve the problem approximately by using 
the equation of motion technique to obtain the exact Green's function g^ a of 
the isolated interacting system, and calculate G^ a from g^ a by using the non- 
interacting tunneling self-energy U of the noninter acting Dyson equation Eq. 
(3.18). This is actually equivalent to performing a Hartree-Fock factorization 
of the higher correlation functions. 

For simplicity, we denote the Green's function of the isolated dot (governed 
only by the Hamiltonian He) by 

g* <T (t) = -id(t)({d a (t) : dU0)}) 

where the subscript (• • -}o indicates that the average is taken for the isolated 
dot. By taking a time derivative, we get 

• =S (t)+e(t)({d*(t),dU0)}} (3.21) 

We use the Heisenberg equation of motion to obtain d a , 

id a = [d a ,H c ] = e a d a + U[d a ,n^di] = e a d a + Ud a n^ (3.22) 
where a = —a. Substituting this result in Eq. (3.21), we get 

= S(t) + e a g^(t) + UgW(t) (3.23) 
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where we define a 2-particle correlation function 

ffSW^W(K(*KW,4(o)})o 

We can also write down its own equation of motion, 
<% (2) 

i^ = S(t)({d a n 9 ,4})o+6(t)({d a (t^ 

Since [n^,d a ] — 0, the first term is simply 6(t)(n^}. The last term is zero 
since n a commutes with He- In the second term we substitute the result of 
Eq. (3.22), so that we finally get 

ft ( 2) 

lJ dT = mM + e ° 9 ™ (t) + Ug ™ {t) (3>24) 

In Fourier space, Eqs. (3.23) and (3.24) become 

(u + i5 - e a - U)g£>(u) = (n a ) 
Substituting the second equation in the first, we get 

(n a ) 1 - (n a ) 



UJ + i5 — Err — U 00 



which is the exact Green's function of the isolated dot. The corresponding 
spectral function has two peaks: one peak at uj — e a with spectral weight 
1 — (ria) corresponding to the probability of the site being occupied by a 
single electron of spin a, and a second peak at uj — e a + U with spectral 
weight (ria) corresponding to double occupancy of the site with electrons of 
opposite spins. 

To obtain the nonequilibrium Green's function G^ a of the connected dot, 
we make the following ansatz: we take the interacting Green's function of 
the isolated dot g^ a as the unperturbed Green's function Gq, and add the 
noninteracting tunneling self-energy Eq. (3.19): 



which gives the following result, 

'- iS - e a - [1 - (n a )]U 



(oj + id - e a ){oj + id - e a - U) - (w + i5 - e a - [1 - (ria)]U)S R (uj) 

(3.25) 

In the simplest case, we neglect the energy dependence of the self-energy, and 
have U R — — i{Ti, + -T^)/2. Note that Eq. (3.25) is actually a self-consistent 
equation, if we take into account the change in spin population (n a ) out of 
equilibrium. Indeed, (n a ) is given out of equilibrium by 
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(n a ) = (dUa) = J |^G<» 

where the lesser Green's function of the dot is given by the Keldysh equation 
which again involves G R a (uj), 

G£» = G; ff ( w )^( w )Gi( w ) 

where U < is the lesser noninteracting tunneling self-energy Eq. (3.20). 

Once G R a {uo) is obtained by a numerical solution of the self-consistent 
equations, we can calculate the current J and the differential conductance 
dJ/dV as a function of bias voltage V — (/ic — HR)/e. We can also add a gate 
potential V g such that the single-particle energies are shifted, e a — > e a + eV 9 . 

In the following figures, we illustrate results for F L — — 0.02 eV, 
e-j- = —0.65 eV, e| = —0.45 eV and £/ = 1 eV at zero temperature. The 
chemical potentials are defined symmetrically fiL = — eVb/2. The figures 
illustrate clearly the Coulomb blockade effect. 

In Fig. 3.1, we display the Coulomb blockade diamond plot, the differential 
conductance dJ/dV as function of gate voltage and bias voltage. Note that 
linear response theory would give only the Vb = line whereas with the 
nonequilibrium theory, we get the full bias-gate voltage diagram. 

In Fig. 3.2, we plot two cuts of Fig. 3.1 along V g , for two values of bias 
voltage: Vb — (equilibrium case for which the differential conductance dJ/dV 
is just the equilibrium two-terminal conductance G), and Vb = 1 V (out of 
equilibrium: the equilibrium conductance peaks split due to the applied bias). 
Finally, in Fig. 3.3 we plot the current J itself as a function of V g and Vb- 



3.2 Linear Response for Steady-State and Homogeneous 
Systems 

We now consider the limit of linear transport, that is, we keep only terms to 
first order in the electric field. We will study the QBE Eq. (2.69) in the 
absence 3 of magnetic field B = 0, and for steady-state (d/dT — 0) and 
homogeneous (Vr = 0) systems. Consider first the equation for the retarded 
Green's function Eq. (2.68). Neglecting the possible dependence of the self- 
energy on the electric field which is a small effect, we have 

G R = - / yR +0{E 2 ) (3.26) 

which to linear order in the electric field is just the equilibrium Green's func- 
tion G R (p, lo). This simplifies drastically the solution of the QBE since all the 
retarded and advanced quantities can be taken in equilibrium. The QBE is 

3 The B / case can also be considered, but then the equilibrium Green's function 
has a dependence on the magnetic field which in general breaks translational 
invariance and makes things more complicated. 
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fu=-0.65 eV, e,=-0.45 eV, U=1 eV, r, =r =0.02 eV, T=0 




y 





Fig. 3.1. Coulomb blockade diamond: differential conductance dJ/dV in units of 
e 2 /h, as function of gate voltage V g and bias voltage 14- Note the doubling of the 
lines due to the level spacing Ae = e; — e-f = 0.2 eV, and the large U = 1 eV spacing 
of the conductance peaks. The equilibrium conductance peaks at 14 = split for 
nonzero applied bias 14 ^ (see also Fig. 3.2). 
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e T =-0.65 eV,£j=-0.45 eV, U=1 eV,r L =r R =0.02 eV, T=0 

0.8 1 1 1 1 1 1 i 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 



V Q (V) 



Fig. 3.2. Differential conductance dJ/dV as function of gate voltage V g , both at 
equilibrium (Vj, = 0) and out of equilibrium (V& ^ 0). 
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-2 -1.5 -1 -0.5 0.5 1 1.5 2 
V Q (V) 



Fig. 3.3. Current J in units of e/h, as function of gate voltage V g and bias voltage 
V b . 
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eE-j y--^^-) • Vp + (v p + VpRei7) 



d_ 
duj 

dU< „ dReG 



VpReG+^-Vp^j = U < A-TG < (3.27) 

This equation can be drastically simplified in the linear response regime since 
the terms in the bracket on the left-hand side can be taken in equilibrium, 
since they are already multiplied by the electric field. This term can then be 
simplified by using the well-known equilibrium expressions, 

Geq(p 5 w ) = in F (w)A(p, uj) , £^(p, uj) = in F (uj)r(p, uj) , 

A(p ' w) = ^TTW ReG{p ^ ) = ^TfW ' (3-28) 

where cr(p, uj) = uj — e p — Re U. In Eq. (3.27), all terms proportional to n F (uj) 
vanish and we are left with 



| A 2 (p, ^eE • [(v p + V p Re E)r + aV p r] = U < A - rG< (3.29) 

UU) 



Using equation (3.29) to derive transport coefficients is equivalent to using 
the Kubo formalism. 

3.2.1 Example: Conductivity from Impurity Scattering 

To illustrate how transport coefficients are calculated from the QBE, we study 
the simple example of the conductivity from impurity scattering in the di- 
lute limit. The left-hand side of Eq. (3.29) contains equilibrium quantities as 
discussed before. The right-hand side (collision term) vanishes identically in 
equilibrium as is immediately seen from the relations (3.28). We keep only 
the terms O(E) in the collision term. As seen in Eq. (3.26), retarded quan- 
tities and the derived quantities A — —2ImG R and r — —2ImS R have no 
O(E) term, hence they are taken in equilibrium. Consequently, only the O(E) 
terms in the nonequilibrium functions U K and G K in the collision term are 
needed. Consider first the collision term. We consider the self-consistent Born 
approximation in which the usual equilibrium time-ordered self-energy is (Fig. 
3.4) 

S(p,uj) = n i Vo + n i J |yp_ q | 2 G(q, uj) (3.30) 

The nonequilibrium contour-ordered self-energy has the same structure, hence 
the nonequilibrium lesser self-energy is 

S<(p,uj)=n i J ^|V p _ q | 2 G<(q,u;) 

Consider expanding the nonequilibrium Green's function G < (p, uj) as the equi- 
librium piece Gf q (p,uj) plus a nonequilibrium piece Gf(p,uj): 
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q = + / 



Fig. 3.4. Self-energy for impurity scattering in the self-consistent Born approxima- 
tion. 



G < (p, uj) = in F (uj)A(p, uj) + G< (p, w) (3.31) 

then the equilibrium part of the collision term vanishes as mentioned earlier, 
and the collision term is 

U<A - rG< = E<A - rGf = I[Gf] 

where 



We now make the following ansatz for the nonequilibrium part Gf, 

G< (p, uj) = -iA(p, uj) ^eE • v p A(p, uj) (3.32) 

where A(p,uj) is a function to be determined. The collision term is then 

drip 



I[G<]=iA(p,uj)r(p,uj) 



duj 



xeE- (v p A(p,uj)- j^J) J (^3l y p-q| 2 ^(q^) v q^(q^)) 

(3.33) 

We now consider the left-hand side of Eq. (3.29). In the dilute limit, we 
neglect terms involving products of r and ReU since they are 0(nf). Hence 
we have 
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(v p + V p Re u)r + aV p r = v p r e p ) v p r + 0{n\) -> v p r 

where in the last step we have neglected the term proportional to u) — e p since 
it is multiplied by A 2 (p,tv) in Eq. (3.29) which is a strongly peaked function 
around 10 — e p in the dilute limit (Re U oc rii is small). The QBE thus becomes 

| A 2 (p, u>) ^eE • v p r(p, w) = iA(p, w)r(p, w) ^ 
xeE- ^VpA(p, W )- -^-y y ^3|y p _ q | 2 A(q,o;)v q A(q, 

from which we extract an integral equation for the unknown function A(p, lv), 

v p A(p,u) = |A(p,o;)vp + * ^ y ^-p|Vp_ q | 2 A(q,w)v q yl(q,£x;) 

which is reminiscent of the integral equation satisfied by the vector vertex 
function in the ladder approximation of the diagrammatic Kubo analysis. In 
fact, the approximations we have made on the QBE are equivalent to the 
ladder approximation in diagrammatic language. 

The electric charge current j is given in general by 

j(R,T) = -ie||^v p / ^G<(p, W ,R,T) 

From the decomposition Eq. (3.31), the equilibrium piece will give a vanishing 
contribution and only Gf in Eq. (3.32) will give a nonvanishing contribution 
to the current. We therefore obtain 

from which the conductivity tensor is easily extracted, 

,2 



e 

m 2 



f d 3 p f dw ( dn F \ ... 



which is equivalent to the Kubo formula result. 



3.3 One-Band Electrons with Spin-Orbit Coupling 

Consider a generic translationally invariant spin-| Hamiltonian with spin- 
orbit coupling, 

C-zV) 2 

H {-iV) = v - ; 1 + A(-zV) • a + /x B B • a 
2m 
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where A(p) is odd in p in general due to time-reversal symmetry, but for 
simplicity we neglect higher-order terms and assume it is only linear in p. 
Consider for example a two-dimensional electron gas with Rashba and Dres- 
selhaus spin-orbit coupling, 



so 



i(Py(J a 



Px& % 



P(Pa 



where a is the Rashba coupling and (3 is the linear Dresselhaus coupling. Then 
we simply have 

A x (p) = ap y + (3p x 

Ay(p) = -(PPx + PPy) 

Orbital magnetic effects are taken care of by the Peierls substitution Eq. 
(2.60). The Green's functions and self-energies in the Wigner representation 
are now 2x2 matrices. The Kadanoff-Baym driving term [Gq 1 — U, G < ] p , ti ;,R,T 
in the presence of electric and magnetic fields can be derived from the general 
expression Eq. (2.62). We obtain 



[G 1 - U, G < ] p ,w,r,t = i 



/ 

+ v D • [ V R + eE— ) + e(E + v p x B) • V p 



-[(A + fi B B) ■ a, G<] + | {V P (A • ^),V R G<} 

dG< 



+ feE-<^ Vp(A-o-), 



duo 



■ [Vp(A-o-) x VpG* 



(3.34) 



where v p = p/m as before, we write simply A = A(p) and we use the notation 
of Eq. (2.65). The anticommutator is derived in a similar manner from Eq. 
(2.63), 



|{^0 1 _ U, G ! ' R }p,a;,R,T — 



c-ep + ^(v R + eEA +e BxVp 



G 



-|{(A + /i fl B) • a, G R } + j[V p (A • <r), V R G fl ] 
dG R 



V p (A-<t) 



eB ■ {V P (A • a) x V P G K } 



(3.35) 



The terms on the right-hand side of the QBE Eq. (2.66) and the equation 
for the retarded Green's function Eq. (2.67) are still given by the general 
Kadanoff-Baym commutators Eq. (2.64). Equations (3.34) and (3.35) are valid 
to arbitrary order in E and B. 

In the case of homogeneous (Vr = 0) transport in zero magnetic field 
B = 0, the retarded Green's function picks up a contribution linear in the 
electric field, as opposed to the case of spinless electrons in section 3.2 where 
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the leading order was 0(E 2 ), Eq. (3.26). This is because the commutator in 
Eq. (3.35) 

" ^ , dG R ' 
Vp(A-o-), 



^0 



does not vanish in general, hence the solution of Eq. (3.35) to linear order in 
E is not just the equilibrium Green's function. 



3.4 Classical Boltzmann Limit 

In this section, we derive the classical Boltzmann equation (2.1) from the QBE 
for 1-band spinless electrons, Eq. (2.69). First, we assume that the quantities 
Re U and Re G are essentially constant so that their derivatives vanish. In the 
QBE, these functions renormalize the transport coefficients that are found 
from the classical Boltzmann equation. We also perform the inverse of the 
change of variables Eq. (2.40), so that we are back with the untransformed 
frequency i?. Hence we have 



+ v p • V R + e(E + v p x B) • V p 



G < {p,f2,R,T) = r > G < - G > U < 



" d_ 
dT 

(3.36) 

Consider Eq. (2.25) for the spectral function. It has been shown by Kadanoff 
and Baym that Eq. (2.25) for the nonequilibrium spectral function can be 
solved in the gradient approximation by the following ansatz, 

A(p,tf,R,T) = 



{Q - e p - ReU - U(R)) 2 + (r/2) 2 

where here U(R) — — eE • R. In the quasiparticle approximation, we assume 
that r is very small so that the spectral function approaches a delta function, 

A(p, n, R, T) = 2tt5({7 -e p -ReU- U(R)) (3.37) 

It is well-known that for equilibrium Green's functions, the fluctuation- 
dissipation theorem holds, 

G<(p, tt) = iA(p, Q)f{Q) (3.38) 
G>(p, Q) = -iA(p, f2)[l - f{Q)} (3.39) 

where f(f2) is the distribution function. Out of equilibrium we define a non- 
equilibrium quantity F(p, f2, R, T) such that the nonequilibrium Green's func- 
tions G <,> (p, Q, R, T) satisfy analogous relations, 



G< (p, Q, R, T) = iA(p, H, R, T)F(p, f2, R, T) (3.40) 
G > (p,/2,R,T) = -iA(p,n,K,T)[l - F(p, i?, R, T)] (3.41) 
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Note that these relations are consistent with the exact relation G K — G > — iA 
and merely constitute a definition of F. However, within the quasiparticle 
approximation Eq. (3.37), the Q dependence of F is seen to be redundant 
because of the delta function and it is sufficient to consider a three-variable 
function /(p, R, T) = F(p, f2, R, T) which plays the role of a nonequilibrium 
distribution function. Equations (3.40,3.41) thus become 

G<(p, Q, R, T) = 27Ti5{f2 - e p - Re U - C7(R))/(p, R, T) (3.42) 
G> (p, n, R, T) = -2m5{Q - e p - Re £ - C/(R))[1 - /(p, R, T)] (3.43) 

where it is seen that the correlation functions G <,> contain information about 
the energy spectrum in addition to information about the nonequilibrium dis- 
tribution. For a general spectral function A(p, f2, R, T), the ansatz (3.40,3.41) 
with F(p, Q, R, T) = /(p,R, T) independent of i? is called the Kadanoff- 
Baym ansatz [8]. We now substituting Eqs. (3.42) and (3.43) in Eq. (3.36) 
and integrate over all i?, which is easy because of the delta functions. We 
obtain 

(^+v p .V R + e(E + v p xB).V p j/ = I[f] 

which is the classical Boltzmann equation (2.1), where the collision term /[/] 
is 

/[/] = -i£<(p, H, R, T) I fi=£p+Re E _ eE , R {l - /(p, R, T)] 
-iZ>(p, n, R, T) | fi=£p+Re £ _ CE . R /(P, R T) 

We can also perform again the Mahan-Hansch transformation Eq. (2.38), 

/[/] = -i£< (p, uj, R, T) | _ £p+Re E [1 - /(p, R, T)] 

-i£> (p, uj, R, T)l =ep+ReE f(p 7 R, T) (3.44) 

which can also be written simply as 

i[f] = s in [i-f]-s out f 

where S in — —iU < is an in-scattering rate and ,S out = zX 1> is an out- scattering 
rate. 

For the sake of definiteness, we investigate a simple example, where colli- 
sions are due to electrons scattering off static impurities in the dilute limit. We 
treat the interaction at the level of the self-consistent Born approximation. At 
this level the self-energy consists of one diagram with one resummed fermion 
line and two connected impurity lines. The equilibrium result for time-ordered 
quantities is well-known (Eq. (3.30)), 

E(p,w) = riiVo + rii \ — — 3|V p _ q | 2 G'(q,a;) 
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where Vk is the Fourier transform of the impurity scattering potential and rn is 
the density of impurities. Since diagrammatic perturbation theory is formally 
the same in or out of equilibrium, the relation still holds for the Keldysh 
contour-ordered functions and the real-time components follow by the usual 
analytic continuation (which is trivial here since there are no products of 
contour-ordered quantities): 

r<'>(p,cu,R,T) = n, J ^3|K p _ q | 2 G < ' > (q,w,R,r) 

Substituting into the collision term Eq. (3.44) and using Eqs. (3.42,3.43), we 
obtain 

I[f] = 2™, J ^|Vp-ql^(e P " e q )[/(q,R,T) - /(p,R,T)] (3.45) 

which is the standard result for impurity scattering. For weak electric fields, 
this collision term can be brought in the following form. Assume a homoge- 
neous system. To linear order in E, the distribution function is 

/(p,T) = /o(p)+C(p,T)p.E (3.46) 

where /o(p) = /eq(e p ) is the equilibrium distribution and C(p,T) is a scalar 
function. Then 

/(q,r)-/(p,r) = c(p,r)(q-p).E 

since the delta function 5(e p — e q ) in Eq. (3.45) enforces e p = e q and p — q. 
Hence the collision term becomes 

I[f] = -C{p,T)2vn % J -0^ |V p _ q | 2 <S(e P - e q )(p - q) ■ E (3.47) 

Choose p along z, and parameterize E with (9, <fi) and q with (9', <fi'): 
p=pz 

E = ^(sin 9 cos <fi x + sin 9 sin (j) y + cos 9 z) 
q = p(sin 9' cos <j)' x + sin 9' sin <fi' y + cos 9' z) 

where we see that cos^' = p • q = cos^ P)q . We get 

(p — q) • E = pE(cos9(l — cos9') — sin^sin^' cos(0 - (j)')) 

In the integral over q, the term cos(0 — (j)') will vanish upon integration over 
(j)' , hence it can be omitted. Inserting the remaining first term in Eq. (3.47) 
and using Eq. (3.46), we get 

/(p,r)-/ (p) 
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which is the usual form of the collision term in the relaxation time approxi- 
mation, and the inverse relaxation time is 

1 f d 3 q 

— = 27m; / 77^-T^|Vp-q| 2 ^(e P - e q )(l - cos# p , q ) 

Tp J (27T)' 3 

Notice the important factor (1 — cos# P;q ) which describes the effectiveness of 
large-angle scattering in destroying momentum. 
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